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ABSTRACT 


Filmwise condensation of steam on externally-finned tubes 
is a very complex process. Recent experiments have shown 
that enhancement ratios (ratio of steam-side heat-transfer 
coefficient to that of a smooth tube having the same diameter) 
exceeded the area enhancement produced by the fins. Moreover, 
the enhancement ratios for fully flooded tubes exceeded the 
values predicted by a simple, one-dimensional conduction 
model by a factor of 2 to 4. A new two-dimensional conduction 
model was developed, which showed that the one-dimensional 
model overpredicted the two-dimensional results for high 
conductivity tube-metals such as copper by as much as 133%. 
The two-dimensional model also showed that variations in fin 
thickness or spacing can result in an overprediction by the 
one-dimensional model of the two-dimensional results by as 


much as twenty-two percent. 
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A. BACKGROUND 

Because the size and weight of main propulsion and 
auxiliary systems are inversely proportional to the ease in 
which heat or energy is exchanged, a major thrust of the 
bescarech and development program of the U.S. Navy is toward 
the design and development of equipment that perform at 
fegaer eriuelency with a reduction in size and weight. In 
particular, one area of great interest is the reduction in 
size and weight of steam condensers. 

Current equipment is of the bare or smooth tube type of 
heat exchanger in which the steam is condensed on the outside 
of the tubes (the shell side). The heat-transfer effective- 
ness is governed by the amount of surface and the overall 
resistance to the flow of heat. This overall resistance, in 
turn, is the sum of the vapor-side, tube wall, and water-side 
thermal resistances. Reduction of any or all of these 
resistances results in an increase in heat-transfer performance 
with an associated reduction in condenser surface area which 
translates into decreased size. 

Both the vapor-side and tube-side thermal resistances can 
be considered as the reciprocal of the so-called hS product 
where h is the coefficient of heat transfer and S is the 
surface area through which the heat flows. Because of the 


curved surface of the exchanger tubing, the vapor-"space" 


1G 


and tube-side surfaces are not equal and an adjustment must 
be made to the tube-side heat-transfer coefficient to allow 
for a design on the basis of the vapor or shell-side surface. 

It is a fact that the vapor-side hS product (reciprocal 
of resistance) can be as high as two or three times the tube- 
side hS product. Thus, it would appear reasonable, in an 
attempt to decrease the overall resistance, to attack and try 
to increase the hS product on the tube side. Unfortunately, 
little can be done in the arena of increasing h significantly 
on either side although the vapor side has greater possibili- 
ties in this regard and, by its very location, not much can 
be done with regard to the increase of the tube-side surface. 

It is for the foregoing reasons that recent efforts have 
focused on the use of externally-finned tubes in steam con- 
densers to increase the flow of heat by increasing the vapor- 
side hS product [Refs. 1,2,3]. 

For the past few decades, it has been accepted that the 
use of externally-finned tubes would be inappropriate mainly 
because of the liguid retention (flooding) that occurs between 
fins (Figure 1.1). This flooding is a result of the relatively 
high surface tension of water (at 50°C water has a surface 
tension of 0.068 N/m, while refrigerants possess values about 
0.015 N/m). One of the unavoidable results of flooding is the 
considerably lower heat transfer through the flooded portion 
than that through the unflooded portion. For the case of an 


externally finned copper tube with steam as the condensate, 
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Tube Wall 


Retained Condensate 


Ficure l.l. Schematic of Condensate Retention on Finned Tubes 
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the copper has a thermal conductivity 600 times greater than 
that of water (385 W/m-K for copper compared with 0.6 W/m-K 
for water). Due to its overall importance in the heat- 
transfer scheme, the flooding of the lower portion of the tube 
by condensate has been studied theoretically and experimentally 
by several investigators [Refs. 4,5,6,7,8]. These studies have 
shown that the flooding angle depends mainly on the fin 
spacing and the ratio of surface tension to density of 
condensate as given by the following equation for the case of 
rectangular-section fins [Ref. 5]: 

Ik 


bo = cos” [1 - 38] (1.1) 





Recent work by Yau et al. [Ref. 9] and Wanniarachchi et al. [Ref. 1] 
has shown that even in the cases with considerable flooding, 
the heat-transfer performance exceeds that expected on the 
basis of area increase due to the fins. 

In 1948, Beatty and Katz [Ref. 10] developed the first model 
to predict the condensation heat-transfer coefficient for 
horizontal, externally-finned tubes. Their model ignores 
the existence of flooding. Although their model agrees to 
within 10 percent for condensing fluids possessing low surface 
tensions, it tends to be rather inaccurate for high-surface- 
tension fluids, such as water. Several other investigators, 
Rudy and Webb [Ref. 2] and Owen et al. [Ref. 8] attempted refinement 


of the Beatty and Katz model but still were unable to predict 
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the heat-transfer coefficient to a reasonable degree of 
accuracy. 

Rudy and Webb used the Beatty and Katz model for the 
unflooded portion, while neglecting heat transfer through the 
flooded portion. With this modification, Rudy and Webb under- 
predicted data for low-surface-tension fluids by up to 
50 percent. Owen et al. attempted to improve the Rudy and 
Webb model by allowing heat transfer through the flooded 
portion. In order to compute heat transfer through the 
flooded portion, they assumed parallel conduction paths 
through the fins and the condensate that is flooded in the 
interfin space as shown in Figure 1.2. With this modification, 
they improved the predictions slightly; underprediction for 
low-surface-tension fluids decreased to about 40 percent. 
However, as shown by Honda and Nozu [Ref. 11],both of these 
models overpredicted the steam data of Yau et al. by a factor 
of up to 3. The fact is that the use of Beatty and Katz 
model, which does not take surface tension forces into consi- 
deration for the unflooded portion, may be responsible for 
this large discrepancy. 

The most recent and accurate attempt at modelling the 
problem was accomplished by Honda and Nozu [Ref. ll]. Their 
besumting solution, a LouUrth—-order differential equation for 
the film thickness, was solved numerically subject to reasonable 
boundary conditions. An agreement to within +20% was achieved 


by their model in predicting data experimentally obtained 


Jig) 
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Figure 1.2. One-Dimensional Parallel Conduction Model 


Ales: 


for 11 fluids and 22 tubes. However, the disagreement with 
steam condensation data was as hich as 40%. 

To understand the condensing heat transfer performance 
more clearlv, it 1S convenient to define the enhancement ratio 
as proposed by Wanniarachchi et al. [Refs.1,2]. The enhancement 
ratio is defined as the ratio of the steam-side heat-transfer 
coefficient of a finned tube to that of a smooth tube with the 
Same diameter as the finned tube root diameter at the same 
heat flux. 

The present investigator has been mainly motivated by the 
heat transfer through a finned tube that is completely flooded 
by condensate. As shown by Winniarachchi et al. [Ref.2], afinned 
tube with a fin spacing of 0.5 mm, and a height and thickness 
of 1.0 mm and 1.0 mm respectively showed a considerable 
improvement in heat-transfer performance even thouch this tube 
was completely flooded by condensate. In fact, this tube 
showed an enhancement ratio of about 2.5 under vacuum (85 mmHg) 
and 4.5 at atmospheric pressure. These ratios even exceeded 
the area enhancement (A;/AL) » which was 2.11. If one- 
dimensional conduction is assumed (see Owen [Ref. 8]), the expected 
enhancement ratio was computed to be about 0.8. This unex- 
pectedly high enhancement ratio may be, at least in part, 
explained by two-dimensional conduction and convection in 


the interfin space. 


fies 


Be 


tb 


OLIRCTive. 


The main objectives of this thesis are: 


Develop a two-dimensional numerical model of a horizontal- 
finned tube fully flooded with steam condensate. The 
program using the model shall be interactive to allow 
different cases of fin thickness, fin spacing, and 
tube-metal thermal conductivity to be studied, 


Using a thermal analyzer computer program, determine 
the temperature distribution plots and heat flux for 
each case, 


Compute enhancement ratios for each one of the cases 
studied and compare these ratios with those computed 
by using a one-dimensional model, 


Compare the one and two-dimensional enhancement ratios 
and comment on the importance of the role of two- 
dimensional heat transfer in the overall heat-transfer 


performance, 


Determine the value and usefulness of the two-dimensional 
model as an aid in the prediction of enhancement ratios 
for various cases. 
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II. DEVELOPMENT OF THE TWO DIMENSIONAL MODEL 


Peeeelae PRINCIPLE OF THE FINDTTE “DIFFERENCE 
1. Establishment of the Spatial Increment 
Consider the three-dimensional Laplace steadv-state 


equation for the spatial temperature distribution in a 





Sent gurdtlLon 
2 2 2 
A aay A) Bak ao a 
= + ; 5 = 0 (2A1.) 


It may be solved analvtically, subject to certain boundarv 
conditions, and the solution is called an analytical solution. 
The finite-difference solution is arrived at numerically and 
applies to a network of points called nodes, with the 
temperatures at the nodes being estimated. The numerical 
solutions obtained by analytical and finite-difference 
methods can be made to agree with any desired accuracy. 
Consider a grid of nodes. Applications of equation 
(2.1) at a particular node whose temperature is TS allows 
estimation of the second derivatives in terms of the 
temperature To and the six related temperatures at the 
adjacent nodes of the grid. In Figure 2.1 an arbitrary 
function T = f(x) is shown which considers only the trace 
of T = £. (x) even thouch eouation (2.1) shows that there 
will also be a functional relationship between T and yv, 


namely, T = f, Cy) and between t and z, namely, T = f(z). 
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°C 
4 ec 


Temperature, 


x Xo Xo 
Distance x, M 


Figure 2.1. Approximation of the First Derivative 


ZZ 


If y and z are held constant, the first derivative 
of equation (2.1) may be approximated by the change in T 


divided by the change in x: 


on = ger Gayl pee Mote op ata ied 
X AX, + Ax. AX, + Ax, AX) + Ax. 
(2.2) 


An estimate of the second derivative is obtained by 
estimating the first derivative in each of the two intervals 
and dividing the difference of the first derivatives by the 
total interval over which the first derivative estimates 
apply. The total interval will be the distance between the 


centers of the two intervals, AX, and AX.- Hence 


5p ; AT,/AX, — AT)/AX, 
de i 1 
oxX > AX. =- 5 AX, 
at AX, = AX. = Ax, 
5 ep (T,-T) /Ax - (LT -T)) ZAx _ Ty - T. - 2T —— 
yx AX ax 
A similar development using T = £, ly) with temperatures 
T, and Ty yields 
5 ep Z 1 ou Ze 
——S SO (24) 
AY 
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and for T = £,(2) with temperatures T. and Te 


5 i 2 =~ = (2.5) 
OZ 


"L* 727 *'o , Bg 4 oe 8 Ue eee 
Ax Ay Az : 


and if Ax = Ay = Az, the result is 


Further justification of the result indicated by 
equation (2.3) comes from a consideration of a Taylor Series 


Expansion. Recall that the Taylor Series Expansion is given 


by 


f (ae Ea) (x-a)° -- eee) ~ om teas 








+ Elica? + oo. + SEB (xa)? (25 


Two cases will be examined. First for f(x) = T(x) 


at x == Xi then 


| 
FF 


T(x=x,) aa 


II 
a 


T(x=xX)) 


24 





2 2 
a at lacie iedeer 3 
Gx ' x ax xX 
Oo O Oo 
4 
A 
Ta aa (2.8) 
24 4 
ax “ls 
Oo 
Bheonmeror f(x) = T(x) at.x = Xo 
T (x=x,) = T. 
T(x=x) = iba 
Qa = 7s 
X-a = -AX 
With these in eauation (2.7) 
- 2 3 
Ty =t,- g| ax + 7 <5 ax” - =5 4) ax 
Ax“! x ax x 
Oo Oo 
4 
eee | oe (2.9) 
24 
fee Xo 


When equations (2.8) and (2.9) are added, one obtains 


1a‘ 


“t 
5 1D —q (Ax) : Seer G20) 


Zo 


If Ax is small, then all terms above the second order 
in Eq. (2.9) may be considered negligible. Thus upon 


rearrangement, it is seen that 


aT T, + T, - 27 
dx ' 7 5 (228) 
O Ax 
which is a result that is identical to eaqudtvon (2.3) oe 


course a Similar development will provide a substantiation 
of equations (2.4) 275) Mand 2-6). 
2. Node to Datum Analysis 

Using the approximations for the second derivative, 
a node analysis can be performed to predict the temperature 
distribution in a mathematical model of a physical configura- 
tion. With a node analysis, the configuration of interest 
is divided into small but finite subvolumes with each sub- 
volume being considered as isothermal. It may be then 
assumed that the geometric center of each subvolume is 
representative of that subvolume and this geometric center 
1s called the node-point or node. 

A node analysis consists of an application of the 
First Law of Thermodynamics to generate energy balances for 


each node. This is done for each node by equating 


HEAT IN + HEAT GENERATED = HEAT OUT + HEAT ACCUMULATED 2-51-29) 


2G 


For the two dimensional case displayed in Figure 2.2, the 
development of the node equations for node point (j,k) is 
Shewn in BPigure 2.32 ° Assuming Heat flow by conductiom and 

a nodal dissipation, q, a steady state or First Law energy 
balance requires that the sum of all heat flows (Watts or 
BTU/hr) leaving the node equals zero. Thus with a heat input 


as a negative heat flow outward, 


oh) i 19) aH 43 ov I Tales = Q (2 4a83)) 


In terms of the temperatures in Figure 2.4 and in terms of the 


Fourier Law 


KA (AT) (2.14) 


4] A 
es Se cies 
q3 = a 6 (TT. 
a = Se te, 
and when these are inserted into equation (2.13), the result 


TUS 


2] 





Two-Dimensional Case 


Figure 2.2. 
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Bigure 2.3. Node Equation Point (j,k) 


Zo 


OT, 


qo 
Ly 
7 
(a. 7 
, ee 
og q3 
| 
* q 
"OTs 


Figure 2.4. Node Equation 
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KAY = KA x é KAYy " 
ee a at eet 


KA X = 
sf ag oer T,) =? Gk. (2-15) 


If two conductances are defined by Kl = KAyé/Ax and 


hee= KAXS6/Ay then 


ee me oo ered (T= De) KR oe qs 
Zo) 
and if there is no heat generation (q.=0), the equation 
becomes 
T, (2K1 + 2K2) - K1(T, ) - K2(T.) - K1(T.) - K2(T,). = 0 
(2 Ey) 


For the case where there is heat flow by both 
conduction and convection aS in Figure 2.5, with a heat 


input as a negative flow outward, 


gear q5 ale q, °° dc. > oa 0 (2178) 


For the conduction terms, Gy do and q3) equation (2.14) is 


still applicable. However, for the convection at the boundary 


as indicated by aie Oput(opulbat=yey = 5))) 


Sy) = BS, (2.19) 


a 





In terms of the temperatures in Figure 2.5 





KCK) : 
_ KA 7 
dp eae © 8G D 
q = K(Ax/2 
3 Ay on Ml T 3) 
° (ae hAy 6 (T,-T.) 
and when these are inserted into equation (2.18), the result 
is 
K (Ax/ 2) 7 KAy : K(Ax/2) : 
Ay 6 (TO T,) + ne S(T) T.) + ie 6 (T. T3) 
+ hAys (T\-T.) = ck (2520)) 
Letting Kl = ae K2 = Ae and K3 = hAys, then 
K1(T-T,) + K2(T-T.) + K1(T -T,) + K,(T,-T.) ace 
C22 1) 
which can be simplified to 
dee ak2 KS te KIT, akon KIT, = qs (2a?) 


SS) 


and once again if there is no heat generation (q, =0), the 


equation becomes 


CZK Ee inc ney K3)T) = Kil K3T., = Kil — K3T. 


a 3 


(2 32 


For the case of a node on an adiabatic boundary with no 


heat generation aS in Figure 2.6 


aT + a5 + qa 0 


(20249 


In terms of the temperature in Figure 2.6 and in terms of the 


Fourier Law 


KUAX/ 2) 


cy. = 
Sy = ee é(T. T.) 

Bein ax/ 2) a 
3 iy omy 


and when these are inserted into equation (2.24), 





1s 
K (Ax/2) = KAy _ Rie) 
aa Sen) fue oy ece) ey ae 
BS) _ KAyé 
Letting aw woe 6 and K2 = Tee 
(2K1 + K2)T, -K1T, -K2T, -K1T, = 0 


It is easy to see that when one proceeds in this 


manner for all nodes, an nxn set of simultaneous algebraic 
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6(T -T3) 


the result 


(2.239 


(2329) 


ak 
ts 
Ay BY a 


aa 






Ss, 





equations will result. These may be represented in matrix 


form as 


[ce (2.27) 


where [K] iS an nxn matrix of coefficients and [T] and [Q] 
are nxl column vectors representing the node temperatures 
and heat inputs respectively. 

The solution of equation (2.27) is handled by the 


computer program TVSSI [Ref. 12:Appendix A]. 


B. MODEL BUILDER TVIN 

In order to provide input to program TVSSI, it was 
necessary to create a model of the configuration to be 
analyzed. The overall model of the condensing region 
incorporated three major features. Because the system to 
be modelled was an externally finned tube with the external 
space flooded by condensate, the azimuthal portion of the 
model was chosen to be 1/21 radians. This is shown in 
Figure 2.7 and because the areas of each node contain a 7 
term, the use of 1/27 simplified the construction of the node 
equations by permitting cancellation of the m terms. The 
second feature required that the maximum number of nodes to 
be used for the model was 300, due to input restrictions 
(the dimension statement) of the main thermal analyzer. The 


final requirement necessitated that the program be written so 
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as to allow the user to input the variables listed in Table I 
interactively. 

The basic model is shown in Figure 2.8. The upper 
horizontal edge represents the condensate/vapor interface 
while the lower horizontal edge represents the inner wall of 
the tube. Heat-transfer coeffients were specified for the 
inner-tube wall and the vapor-condensate interface, The 
model assumes that the heat is transferred across the retained 
condensate by conduction only. Due to symmetry, the dotted- 
dashed lines at the two vertical edges of the figure represent 
an isolation of a fin-condensate-tube entity into a repeating 
section. These lines, of course, represent adiabatic 
boundaries. Thus, the model assumes that the only heat 
transferred into or out of the model occurs at either the 
condensate/vapor interface or the inner tube wall. 

The first model constructed consisted of 299 nodes with 
13 columns and 23 rows. The outermost rows and columns 
consisted of half node subvolumes. This facilitated the 
placement of the node centers at both the inner tube wall 
and upper condensate/vapor interface in addition to the 
boundary of the outer tube wall and fin perimeter. Nine rows 
were placed in the condensate region above the fin and in the 
region including the fin and adjacent condensate. The 
remaining four rows were placed in the tube wall region. 

In addition the fin region contained 4 columns of nodes 
while the remaining eight columns were placed in the inter-fin 


Space. 
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ENE 


TABLE I 


VARIABLES FOR MODEL BUILDER TVIN 


Temperature of fluid inside tube 


Temperature outside condensate layer 


Inner radius of tube 


Outer radius of tube 


Heat transfer coef. at inside tube wall 


Heat transfer coef. at condensate layer 


Thermal conductivity at tube and fin 


Thermal conductivity of condensing liquid 


Fin height 


Fin width 


Spacing between fins 


Height of condensate above tip of fin 
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The location of these columns and rows determined the 
height and width increments of the various subvolumes, 
depending on their location in the model. For example, if 
the tube wall thickness is four millimeters, and the tube 
wall area contains four rows, the height increment for a 
node in this region would be 4mm/4 rows =1mm. The width 
increments for a node would be computed in identical fashion 
by dividing the size of the particular region by the number 
of columns occupying that region. Figure 2.9 identifies all 
of the different height and width increments utilized by the 


model. 


C. EVALUATION OF RESULTS 
Pehl LCatlon Of Accuracy 
Verification of the model consisted of using the 
temperatures generated by TVSSI (Using TVIN as input) and 
the conductance values employed in TVIN to compute the heat 
flow across three boundaries in the model: 
(1) condensate upper surface 
(2) fin perimeter and tube outer wall 
(3) tube inner wall 
The foregoing heat flows were computed in each case 
by using the node temperatures generated by TVSSI with the 
associated internode conductance values used in the input 


generated by TVIN. These heat flows, of course, must be 
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ne 











equal and it was observed at the outset that 13x23 model 
yielded certain unacceptable discrepancies in the three values. 

Because the 13x23, 299 node model was found to be 
inaccurate, it was determined that, in order to achieve the 
requisite accuracy, slight adjustment of the number of rows 
and columns was required from run-to-run. 

TVIN was reprogramed to allow flexibility in the 
number of rows and columns utilized in the model as well as 
placement of the nodes in the fin, tube or liquid regions. 
The first change instituted was to fix the number of rows 
contained in the condensate layer above the fin tip at two 
rows. This action was acceptable because the original height 
of this layer was only 0.08 mm compared with a 1.0 mm fin 
height and 3.175 mm tube wall thickness. Fixing the number 
of rows at two permitted the use of additional rows at the 
fin and tube levels. Program TVIN was changed to inter- 
actively accept a specified number of columns (see Appendix 
A). The number of rows was obtained by making the product 
of rows and columns as close as possible to 300. 

A description of the menu-driven procedure for the 
peleceion and alllocation Gs news and columns for a 296=nede 
model example now follows: 

If the response to the query 'How many columns are 
desired?' is 8, the program will provide 300/8 = 37.5 or 37 
rows. This yields an 8x37 = 296 node model. 

Suppose the response to the query ‘How many rows are 


desired for the tube wall?' is 5. Then the program will 
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allocate 37-5 = 32 rows for the fin and condensate layer 
above the fin. 

However, the number of rows in the condensate layer 
above the fin tip is fixed at 2. Thus the number of rows in 
the fin becomes 32-2 = 30. The result is an 8x37, 296 node model. 
A display of the foregoing model set-up is shown in 
Figure 2.10. 

2. Model Verifier Ql 

As mentioned in the previous section, once an input 
file had been created by TVIN, and the node temperatures 
were calculated by the thermal analyzer, TVSSI, the solution 
was verified for correctness and accuracy. This check 
consisted merely of a comparison of the heat flow (watts) 
across three surfaces in the model. An error of 3.5% was 
established as the maximum allowable for a solution to be 
considered acceptable. The model verifier Program Ql (see 
Appendix B), is an offshoot of the model builder TVIN as it 
interactively accepts the same input data required to build 
the model. Once all input parameters are set, the program 
generates and stores each inter-connecting node conductance 
value. 

In the next step, all node temperatures computed by 
the thermal analyzer are read and stored, by location in the 
model. The node numbers of the condensate outer layer, fin 
perimeter and outer tube wall and the inner tube surface are 


determined for each model via an algorithm which requires as 
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input, the number of columns and rows being used in the model 
along with their placement. In addition to determining the 
nodes on these surfaces, the algorithm also computes the 

node number of each of the adjoining nodes in the direction 

of the heat flow. Computation of the heat flow by program Ql 
is then simply a matter of taking each surface node conductance 
and multiplying it by the temperature difference between the 
surface node and its adjoining node. The summation of all 

the surface node values results in the total heat flow across 


the particular surface. 


D. THERMAL ANALYZER TVSSI 

The computer program TVSSI used in obtaining the temperature 
distributions was adapted from program TVSS2 listed by Kern 
and Kraus [Ref. 13:Appendix 2a]. This adaptation consisted 
of modifying TVSS2 to compute solutions in the SI system and 
accept an input file specifically designed by model builder 
TVIN for each externally finned tube case. This program also 
had to be adapted to the IBM 3033 AP system. 

The program utilizes the Cholesky decomposition scheme, 
and because of the linearization of the radiation terms (not 
used in this study), the program is iterative. A more indepth 
description of the Cholesky Scheme pertaining to the program 
1s discussed by O'Hare [Ref. 12:pp. 48-52]. 

For all the cases which were studied, the following 


parameters were set for the analyzer: 
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1. The initial temperatures at which all nodes were set 
and at which the computer began the iterative process 
was 100°C. 


2. An eventual accuracy of .05 between the final and next 
to last iterations was used. 


3. A radiation coefficient convergence factor Of 0.66667 
between iterations was used. This is required by the 
program even if radiation is not considered as a mode 
of heat transfer in the model. 

4. The maximum number of iterations that the computer was 
allowed to perform was set at 12 to avoid excessive 
computer running time. 

5. A damping factor of 0.8 was set as an initial damping 
factor based on the number of non-linear terms in all 
of the node equation. 

The ability of TVSSI to correctly solve a numerical model 
was verified by providing as input a data set listed by Kern 
and Kraus [Ref. 13:pp. 326-336] for an extended surface 
problem where the solutions were derived analytically. The 
computer generated solution agreed completely with the 


analytic solution thereby assuring confidence of analyzer 


generated temperatures. 


47 


III. RESULTS AND DESCUSS EON 


A. INTRODUCTION 

During this thesis effort, 15 different cases were studied 
using the programs TVIN and TVSSI. Table II lists the param- 
eters that remained constant in all cases. A summary of all 


cases examined and their results are listed in Table III. 


B. SMOOTH-TUBE EVALUATION 

As a method for comparison of the variation in heat- 
transfer performance, the heat flux for a smooth tube with 
the same diameter as the finned-tube root diameter was estab- 
lished as a baseline value. The ratio of the one- and two- 
dimensional heat fluxes to the heat flux of the finned tube 
with fin width approaching zero were designated on and 9. 
respectively. The equation used to determine the finned-tube 


of zero fin width heat flow is given by: 


GO aes) 
Onne O l (3.1) 
iB : D, (in(D_/D,)) ae 
(hj) (Dy) 2s Kh, 


C. ONE DIMENSIONAL HEAT FLOW 

The one-dimensional model utilized in this study assumed, 
as did Owen et al. [Ref. 8], parallel paths for heat flow across 
the tube wall, fin and retained condensate. Furthermore, a 
Simplifying assumption was made that heat flow is purely radial 


and that the heat is transferred across the retained condensate 
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TABLE IT 


CONSTANT PARAMETERS FOR MODEL BUILDING PROGRAM TVIN 


A; = 30000 W/m:K 
Ay = 100000 W/m. kK 
Yr. = §¢6.295Fmm 

a 
ic = Om 2 5) il 

Oo 
& = jie0enun 
Ky = 0.6 W/m-K 
T Se BO aacG 

re) 

= 0 
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by conduction only. This was accomplished by modifying pro- 
gram TVIN in such a way so that the data file it created con- 
tained only those internode conductance values representing 
heat flow in the radial direction while deleting those conduc- 
wamec terms pertaining to axial” heat flow. Programs TVSSI amd 
Ql were utilized to compute the node temperatures and heat 
flow for each case. Once again, a maximum error of 3.5% was 


used when heat flows generated by program Ql were verified. 


D. TWO DIMENSIONAL HEAT FLOW 
The heat flow for the two-dimensional case was evaluated by 


the methods discussed in Chapters I and II. 


E. EFFECTS OF FIN SPACING AND FIN WIDTH ON PERFORMANCE 

Heat fluxes were computed for five different cases of fin 
width and spacing. Three cases had a fin width of 1 mm with 
Spacings of 0.5, 1.0 and 1.5 mm respectively. The other two 
cases held the spacing constant at 0.5 mm, while the fin width 
was adjusted to 0.5 and 0.25 mm. 

1. Effects of Fin Spacing 

This section presents results showing the variations 

bpamieat flux having fin spacing as a parameter. Figures®3-.1, 
3.2 and 3.3 present results for tubes having a fin width of 1.0 
Tieemehcin spacings Gf 0.5, 1-0) and 1.5 mm. “Bach £vqume™repre- 
sents the same spacing variation with the tube-metal thermal 
conductivity being varied from 385.0 W/m.K for copper, through 


45.0 W/m K for copper-nickel, down to 16.0 W/m-K for the case of 
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Secimtes's’ steel. Similarly, Figures 32.4, 3.5uand 3.6 present 
the ratio of the two-dimensional heat flux to the one-dimensional 
heat flux (Q21). For each case of tube-metal thermal conduc- 
tivity, the fin spacing of 0.5 mm shows the best heat-transfer 
performance, while the worst performance occurs at a Spacumg of 
1.5 mm. Figures 3.7, 3.8 and 3.9 present isotherm plots for 
ame “emree cases of thermal conductivity for the £imned ube 
have a 0.5 mm fin spacing and 1.0 mm fin width. Moreover, the 
ratio Q21]1 has its smallest value at a find spacing of 0.5 mm 
and increases with increasing fin spacing. Enhancement ratios 
for the various cases had their maximum value at the 0.5 mm 

fin spacing and this ratio showed a marked decreaase as the fin 
spacing was increased. 

The results indicate several important facts. The first 
and most important is the significant role the two-dimensional 
heat flow plays in this type of extended surface analysis. The 
increasing Q21 ratio with increasing fin spacing shows that for 
the cases of copper and copper-nickel tube-metal, the axial 
heat flow into the fin upper side surface tends to raise the 
temperature of the fin tip resulting in a choking effect of the 
heat flow into the fin tip surface. This is of enormous conse- 
quence because in each and every case, the majority of the 
heat absorbed by the finned tube is accomplished at the fin 
tip. Because the one-dimensional case allows no axial heat 
flow into the fin side surface, a greater temperature gradient 


Occurs at the fin tip reswlermenin a greater heat flow. As the 
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Figure 3./7. Isotherm Plot For Finned Copper Tine 
(Fin and Condensate Regions) 
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Figure 3.8. Isotherm Plot For Finned Copper-Nickel Tube 
(Fin and Condensate Regions) 


60 















e 
. 
a 


sm — =\\ = 
“ \  « acing 


‘ —— bling 
8 — 
ae 5 * 






0.48 6.90 0.8 0.60 0.83 0.70 0.73 


Fiogure 3.9. Isotherm Plot For Finned Stainless Steel Tube 


fin spacing 1S increased, the effect of the fin side surface 

is less pronounced since the poor conductivity of the condensate 
inhibits axial heat flow in the region near the fin space 
centerline. This results in a more cloSer agreement between 

the one-dimensional model and the two-dimensional version. 

This conclusion can be further justified by comparing 
the Q21 ratios for each tube-metal conductivity. The two- 
dimensional effect at the fin tip is decreased as the thermal 
conductivity of the fin material decreases. This does not 
mean, however, that there is no two-dimensional effect for the 
cases with low tube-metal conductivity. As in the case of 
copper, the two-dimensional effect tends to inhibit the heat 
flow into the fin. For the cases of low metal thermal conduc- 
tivity, there is a less pronounced two-dimensional effect at 
the fin tip, but there is a significant two-dimensional heat 
flow into the lower portion of the fin side surface near the 
tube wall. For this reason, the one-dimensional model tends 
to underpredict the two-dimensional model. 

As seen in all the cases for the copper-nickel and 
stainless-steel tube-metal (except the copper-nickel case 
with 0.5 mm fin spacing), all Q21 ratios are greater than 
unity. For the copper-nickel case of 0.5 mm spacing, the reduced 
spacing even with a lower metal conductivity of 45 W/m-K, tends 
to produce a two-dimensional effect at the upper fin tip. 

2. Effects of Fin Width 
The results of the variation of heat flux with fin 


width as a parameter indicate that fin width plays an important 
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role in the overall heat-transfer mechanism. The two cases 

for each tube-metal thermal conductivity investigated were fin 
widths of 0.5 and 0.25 mm respectively with fin spacing held 
constant at 0.5 mm. For each case of tube~metal thermal con- 
ductivity, the fin width of 0.5 mm displayed one of the highest 
heat fluxes, being second only to the case of the 1.0 mm fin 
waren wath 0.5 mm fin spacing. The only significant two= 
dimensional fin tip effect can only be observed in the case 

of the copper tube-metal with a fin width of 0.5 mm. 

The same reasoning as in the previous section can be 
applied in this case. Because the fin tip width has been 
reduced, the total amount of heat transferred into the tip 
region is also reduced and the two-dimensional fin tip choking 
effect has less influence. This concurs with the computed Q21 
rarer Or 1.042" By reducing the fin width further in eacheease 
to 0.25 mm, each of the resultant heat fluxes is less than the 
cases with 0.5 mm fin width. For each case of the 0.25 mm fin 
Teadene che O21 ratio appreaehnedm@l.s, signifying the comparatively 


large two-dimensional heat flow into the lower fin side surface. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A.” “CONCLUSIONS 

1. The two-dimensional numerical model of a horizontal 
finned tube fully flooded with steam condensate developed is 
a useful tool in the study of the relationship between one- 
and two-dimensional heat-transfer effects. 

2. Two-dimensional heat transfer effects play a major 
role in the overall heat transfer mechanism in the case of an 
externally finned tube fully flooded with steam condensate. 

It has been determined that for cases of high tube-metal con- 
ductivity, such as the case of copper, the one-dimensional 
model can overpredict the results of the two-dimensional 

model by as much as 13%, while variation in the fin width or 
spacing can result in the one-dimensional model underpredicting 
the two-dimensional result by as much as twenty-two percent. 

3. Fin width and fin spacing control the magnitude of the 
two-dimensional effect present and the resulting heat-transfer 
performance. An increasing fin width to fin spacing ratio 
will correspond to an increase in the heat-transfer performance. 

4. The one-dimensional model in cases of high tube-metal 
thermal conductivity exhibiting a fin width to spacing ratio 
of greater than unity will overpredict the heat flux predictions 
of the two-dimensional model. 

5. Decreasing the tube-metal thermal conductivity results 


in a degraded heat-transfer performance. 
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6. The two-dimensional fin-tip choking effect is minimal 
in low thermal conductivity tube-metals while the two- 
dimensional heat transfer into the fin lower side surface is 


far more pronounced. 


B. RECOMMENDATIONS 

1. Use the two-dimensional model computed enhancement ratios 
for the various cases to estimate a steam-side heat-transfer 
coefficient and compare with those methods currently being used. 

2. Using the model builderiTVIN and thermal analyzer 
TVSSI, conduct a more in-depth study of the effects of variation 
ef fin width, height and spacing on tne overall heat-transfer 
performance. 

3. Revise the two-dimensional model to include convection 


effects in the retained condensate region. 
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